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We give an explicitly gauge-invariant canonical analysis of linearized quadratic gravity 
theories in three dimensions for both flat and de Sitter backgrounds. In flat backgrounds, 
we also study the effects of the gravitational Chern-Simons term, include the sources, and 
compute the weak field limit as well as scattering between spinning massive particles. 
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Recently, Bergshoeff et al. [l[ found that, in three dimensions, among the class of higher- 
derivative theories defined by the Lagrangian k~ x R + aR 2 + pR^, a special case 8a + 3/3 = 
and ft -1 < (let us call it BHT gravity) and its parity- violating extension, with a gravitational 
Chern-Simons term, have massive ghost-free spin-2 particles in their free spectrum around both 
flat and (anti)-de Sitter [(a)dS] spacetimes. Perhaps, the most interesting feature of the BHT 
model is that it is the first and (apart from some bimetric theories) the only known example of a 
(parity-invariant) theory that provides a nonlinear extension to the Pauli-Fierz mass term for spin- 
2 particles. In addition, being a three-dimensional theory, it is powercounting superrenormalizable 
whose four-dimensional cousin is renormalizable [2J. Therefore, it is possible that the BHT model 
may turn out to be a perturbatively well-defined quantum gravity in three dimensions. But of 
course, unitarity of the model beyond tree level is yet to be checked. 

Various aspects of the theory such as its ghost-freedom and tree level unitarity [H, 0-0] and 
Newtonian limits [5[ have been explored. Also, classical solutions and related issues were studied 



CD m , and supergravity extensions were given in 

In this paper, we give an explicitly gauge-invariant, detailed analysis of the canonical structure 
of the generic quadratic models in 2 + 1 dimensions for both flat and de Sitter (dS) backgrounds. In 
flat space, we also include the gravitational Chern-Simons term in our analysis. It is interesting to 
^ ■ see how at the linearized level BHT theory is singled out as a unique regular "harmonic oscillator" 

(massive free field), which avoids the infamous Ostragradskian instability that ruins every higher- 
time derivative theory [l2T ] . [It was claimed that adding interactions might yield stable higher-time 
derivative theories [13j.] All the other quadratic theories are ghost-ridden higher-derivative Pais- 
Uhlenbeck oscillators at the linearized level. In addition, we discuss the Newtonian limits, 
weak fields, and the tree level scattering of particles with mass and spin in these models. 

The layout of the paper is as follows: Section |TT] is devoted to flat spacetime analysis which 
includes the canonical structure of both the parity-invariant and parity- violating quadratic gravity, 
in addition to the effects of static sources and weak field solutions with circular symmetry. In 
Section IIIIl canonical structure analysis is extended to de Sitter space. Some of the computations 
are relegated to the Appendices. Tree level scattering amplitude between spinning massive particles 
is given in Appendix A. In Appendix B, generic quadratic action is written in terms of two auxiliary 
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fields. Finally, we list some results which may be helpful in the analysis of field equations. 



II. HIGHER-DERIVATIVE SPIN-2 IN FLAT SPACETIME 



We start our analysis of the higher-derivative spin-2 fields in flat space, which is considerably 
simpler than the de Sitter background, which we deal with in the next section. The action 



I 



-H { -R + aR 2 + (3R 2 



(1) 



gives the desired spin-2 model when expanded as g uu = r] uu + h uu , where rj uu is the usual flat 
spacetime metric with mostly plus signature. [Actually, spin-2 here is a misnomer. It should be 
symmetric rank-2 tensor, since without any constraints in addition to spin-2, it has spin-1 and 
spin-0 components. But, in what follows, we will call h au a spin-2 field.] Below, we will also add 
the parity-violating gravitational Chern-Simons term to this action. In practice, to actually get 
the action for huu, it is somewhat more convenient to linearize the full nonlinear field equations 
and then integrate them (after carefully taking care of the overall sign, which will be relevant for 
the discussion of ghosts). Then, the action ([T]) up to boundary terms becomes 



I 



~ / (f'.rh 



Igp + ( 2a + 0) {uTu - d»d v ) R L + pUQ^ 

Kj 



(2) 



Here, the linearized Einstein and Ricci tensors, and curvature scalar read 



\rf v R L , R L = d a dsh a ? - ah, 



- (a^T + d a d u h^ - nh» v - d»d v h) , h = rfK^ 



(3) 



where □ = d a d^ 



-8q + V 2 . Raising and lowering operations are carried out with r] uu . To explore 
the canonical structure and identify the free fields, hav can be decomposed in terms of six a priori 
free functions of (t, x) : 



hij 
hoi 



+ didj) 4> - didjx + Uikdkdj + e jk d k di) £ 



eijdjrj + diN L , 



oo 



N, 



(4) 



where di = di/^J — V 2 . From these, one can compute in terms of three functions: 



'00 



2 ( Cik ® k 



a + di 



(Sij + didj) q - didjcf) - (e ik d k dj + e jk d k d^J a 



where (j) = dcft/dt, etc. Here, q, a, and eft are invariant under gauge transformations S^h au 
d^Qu + d u ( u and are defined as 



q = V 2 N - 2V 2 N L + x, o- = £-V 2 rj. 



(5) 



Note also that eft is gauge invariant unlike the other components of hav ■ Linearized scalar curvature 
is computed to be 



R L = q~Oc 
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Therefore, as required by the Bianchi identity, d^Q't = 0, the number of arbitrary functions reduces 
from six to three. One can use either cj), a, q; or (j), a, Rl combinations. The Einstein-Hilbert part 
of the action can be computed as 

lEH = ~hJ d * x hfiuG ^ = hf d3x + a2 ) ' 

which clearly shows that there is no propagating degree of freedom in the pure Einstein theory. 
To compute the quadratic part, its better to use the self-adjointness of the involved operators to 
rewrite the action as explicitly gauge invariant not just gauge invariant up to a boundary term, 
which will simplify the computations in a great deal: 

W = -^f^ / d 3 x V «□ - VP) R L = J d'xRl 



d 3 xh 



d 3 x 



d 3 x 



/□</> + aUa) 



In the Ip action, the second equality follows after one moves the □ term to and then uses 
and the Bianchi identity. Collecting all the terms, the total action in terms of the gauge-invariant 
combinations is 



I 



1 



d 3 x 



1 



<+(2a + /3) (q - nfi) 2 + pqU<f> 



+ 



d 3 x 



1 2 

up 



(6) 



a describes a single scalar field with mass m 2 = which is nontachyonic for k(3 < and a 
nonghost for (3 > 0, therefore k < 0. For the <ft and q part of the action, the discussion bifurcates 
whether 2a + f3 = 0, or not. Let us first consider the 2a + (3 ^ case, for which the nondynamical 
field q can be eliminated, yielding the action 



Li: 



1 



d 3 x 



13 (8a + 3/3) 
4{2a + P) 



+ 



(4a + ^) 
2k(2o + /3)' 



4k 2 (2a + P) 



(7) 



There are apparently several special points one of which is the BHT limit 8a + 3/3 = 0, for which 
the higher-derivative term disappears. [The 4a + /3 = theory seems special, but it has a tachyonic 
excitation; on the other hand, the j3 = model is ghost and tachyon free for k > 0.] Therefore, 
at the linearized level, the BHT model is actually not a higher-derivative theory, so it escapes the 
Ostragradski instability. The (f) field part of the BHT action reads 



1 

'2k 



d 3 x 



!>□</> + 



1 



which again describes a single degree of freedom with the same mass as a. This is to be expected 
in this parity-invariant theory, since a and <f> are two helicity degrees of freedom of the massive 
spin two field in three dimensions. Also, observe that for cj) to be a nonghost, n has to be negative. 

For generic a and /3, except for 2a + (3 ^ 0, (0 describes a higher-derivative Pais-Uhlenbeck 
14| oscillator which can be rewritten in terms of simple oscillators in the following way. Defining 
new fields as 



9?2 



mi 
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([7]) becomes 



1 



64k(2« + /3) z 



with m 9 given as above and m s as 



(8) 



1 



k (8a + 3^9) 



For 8a + 3/3 < 0, 991 is nontachyonic just like (f2, but unlike it describes a ghostlike excitation 
since its kinetic energy comes with the wrong sign. 



A. 2a 



theory 



We have seen in the above discussion that the 2a + j3 = case is a somewhat singular theory. 
If one naively takes the e = 2a + /3 — > limit in (jSJ, one gets 



1 

8^e 



era: 



/3 



mi 



(□ - mfj <j) 2 - 4e0 (□ - mj) + O (e 2 ) 



which is a degenerate (equal mass) Pais-Uhlenbeck oscillator after a divergent rescaling of eft. But, 
more properly, suppose from the onset at the level of the action, we set 2a + f3 = to get (apart 
from the decoupled a field) 



d?x (^qO(f> — m 2 qcj) 



Variation with respect to <p gives a massive wave equation for q, and vice versa. But, these equations 
do not reveal the ghost structure of the theory. So, let us define q = m 2 g (^1 + ^2), 4> = ^1 — 
which turns the action to 

m 2 J 



I 



dfx 



m 2 g ^f\ - (* 2 D^2 - m 2 g ^l 



Since f3 > 0, ^2 is a ghost excitation. The Newtonian limit of this theory is quite interesting: From 
the general tree-level scattering amplitude computation given in p|, one sees that as in the pure 
Einstein-Hilbert theory, the 2a + j3 = case has a vanishing Newtonian potential between static 
sources: the spin-0 ghost excitation gives a repulsive component which cancels the attractive one 
coming from the spin-2 part. 



B. Adding static sources 

Up to now, we have studied the free field spectrum of higher-derivative gravity. Let us remedy 
this by adding matter with the usual gravity-matter coupling: 

-^source = ~Z j d X h^i/T^ . 



In the case of a static source, T 00 = p(x), T° l = 0, T l i = 0, (in a related context, we somewhat 
generalize this in Appendix A), l sourC e becomes 

Wee = \ J d 3 xNp(x) = ^Jd 3 x f^q + 2A L - ^x) pip) , 
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where in the second equality, we have used the definition of q in ([5]) . After dropping the boundary 
terms and using the symmetry of the Green's function, we have 



1 



1 



- 2 Jdx q — 2 p. 

Redefining ip = <fi + i^^iP and q = q + np, the total action reduces to 



d 3 x 



— (pq — npp + a 2 ^j + (2a + j3) (q — □y?) 2 + (3 (q^p — KpOp — nqp + k 2 p 2 + crder 



Specifically, for 8a + 3/3 = 0, integrating out g, one ends up with 



1 



d A x 



aOa — m 2 a 2 



The last term is the interaction part which gives the attractive (for n < 0) potential energy 



U 



d 2 x p\ 



y z — mi 8tt 



0) 



where we took point sources, pi(x) = mii"(f-^d, P2{x) = 17126^ (x — x^), and -K"o is the 
modified Bessel function. This result matches that of [5|. 



C. Weak field approximation 

It is also highly instructive to capture some of the above results from the nonlinear theory ([1} . 
But, even in the circularly symmetric case, nontrivial exact solutions for which 500 9 rr are n °t 
known, and we have not been able to find one. Nevertheless, since we just need the weak field 
approximation, we can do the following: The ansatz 

^2 = _ f (r) dt 2 + h lkl dr 2 + r 2^ 2j 

can be inserted into the action ([T]), which is to be varied with respect / (r) and b (r) [See the details 
of this Weyl trick in [15j] . For the sake of simplicity, let us just consider the BHT theory. Then, 
an approximate solution can be found by setting /(r) = 1 + f r dra(r), b(r) = 1 + f 7 drv(r), 
where a and v are small. At the first order, we have 

-v + 2/3v" + 2/3a" + r/3a'" = 0, (10) 

K 

2 

f3r 2 a" + -r 2 a + 2r(3v ' - 2f3v = 0. (11) 

Here, ' denotes differentiation with respect to r. v can be determined as v = a + |a'. Putting it 
back to (fTTjl gives 

r 2 a" + ra! -a (m 2 r 2 + l) = 0, (12) 

which is solved by a (r) = c\ I\ (m g r) + c^Kx (m g r). Recall that goo ~ — 1 — f dra(r), and 
g rr ~ 1 + f r dr (2v (r) — a (r)). Thus, for decaying fields c\ vanishes, and the metric components 
become 

500 ~ -1 + cK (m g r) , g rr « 1 + dK\ (m g r) , 

where c and d are constants related to the mass of the source. This is consistent with our earlier 
result ©. 
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D. Higher-derivative gravity plus a Chern-Simons term 

We will now extend the preceding discussion in flat space by adding a gravitational Chern- 
Simons term 16] 



I 



d 3 x \J —g 



(13) 



where eoi2 = 1 ; an d A* i s the Chern-Simons coupling with an arbitrary sig n. [ Without the a, (3 
terms, but with a Pauli-Fierz mass term, canonical analysis was carried out in 17, 1||] Linearization 
of the Chern-Simons part yields 



l cs 



— J d 3 xe m& Qfd^ v = ±-Jd 3 xa(q + 



□c 



The total action in terms of the gauge-invariant combinations becomes 



/ \ I <r\r 



1 



+ a 2 + (2a + P) (q - □(/>) + P (gD0 + <tD<t) + -a (g + □< 



Assuming that 2a + f3 ^ 0, q can be eliminated to yield the action 



1 

+- 

K 



a a a + I — 



1 



rc/3 V/3(2a + /3) 



+ 



1 (4a + 0) 



/fc (8a + 3/3) (D# + j^±4 w 



4(2a + /3) 



2 (2a + P) 



H 2 / u(2a + /3) 

—7^ 

4K(2a + /3y _ 



1 



2/c//(2a + /3) 



For generic a, ft one can diagonalize this action, but it is rather cumbersome and not particularly 
illuminating, so we just consider the 8a + 3/3 = case, 



P 

Ibht-cs = - I dfx 



aUo - [mt + 



1 



/i 2 /3 2 



a 



+ + m 2 , (</>□</> - m 2 ^ 2 ) L . 



Pfi 



To decouple the a, (ft fields, one possible route is to take the Fourier transform of the fields, put 
the Lagrangian in a matrix form, and then diagonalize the matrix. This procedure yields 

Ibht~cs = ~ f d 3 x f^ + D^+ - mi^i + <!'_□<!'_ - m 2 <I» 2 

where the masses read 



2 2 , 

m ± = m + 



± ml + 



g 2/i 2 /? 2 fiP\ 9 V/3 2: 



and the new fields are defined as 



'n + 


(m. 2 ^ — m 2 


7V + " 


AL 


(m 2 . — m 2 





m q <j) ) ' 



\ 



1 + 



»P 



- 2 

(m\ — m 2 J 



m± agree with those of [H, EH. As the +2 and —2 helicity modes have different masses, it is a 
parity- violating theory as expected. In the P — > limit, which is the topologically massive gravity 
with a single degree of freedom [la ]. m+ diverges and drops out, m_ = — \/j,\ / k. 
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III. HIGHER-DERIVATIVE SPIN-2 IN A DE SITTER BACKGROUND 



Now, we will study the canonical structure of the higher-derivative theory in an (anti)-de Sitter 
background defined by the action 



dTxyJ^g 



1 



(R - 2A ) + aR 2 + (3R 



whose linearization about an (a)dS background yields 



d 3 x y/—g h 



aQ% + (2a + P) (gTU - VV + Rl + P - ±gT ' R L 



where a = i + |f a + and \ji 2 is the cosmological constant which is related to a, /3, k and 
the bare cosmological constant A of the full theory as -p = w^qrg] 1 ± \A — 8kA (3a + (5) 
fl^ |. For the sake of simplicity, we will consider the background to be a de Sitter spacetime, but 
since our results will be analytic in £, in the final expressions one can take t — > i£ to obtain the 
results in anti-de Sitter spacetime. [To keep the signature intact, one also needs to Wick rotate a 
space coordinate]. For dS, we take the metric, g^ u , with which all the covariant derivatives and 
raising-lowering operations should be made, to be in the 



ds 2 = [-dt 2 + dx 2 + dy 2 ^j , 

and define the perturbation as 

- £ 1 h 

Linearized forms of Einstein and Ricci tensors, and Ricci scalar are given as 

R] 



R 



- (V CT V M /w + V CT V, V " D V " V»Vvh) , R L = V a V^r -Uh- -gh, (14) 



where □ = V^V^ = jzij^V ^V u . Decomposition of h^ v into "spatial" tensor hij, "spatial" vector 
hoi, and "scalar" hoo is 



f2 



hoi 



t 2 

P 
t 2 

t 2 



[Sij + ViVj 
(Sij + ViVj 



iV, 



where V.; = V^/y— and the covariant derivative is for two-dimensional space with metric 
Hi = fr^ij- Since the two-dimensional space is flat, then V» — > di and di = di/J—d^. Zik is 
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the Levi-Civita tensor for two-dimensional space, which is related with the corresponding tensor 
density e ik by 



hk = \fltik 



tik 



t 2 



tik- 



The convention for is e\2 = 1 (the convention for Levi-Civita tensor density for the upper indices 
is e 12 = 1 naturally with the induced metric). As a result, the final form of the decomposition is 



h t 



Sij + 



did^j 4> - didjx + (eikdkdj + e jk d k d^ £ 



I 2 i 2 
hoi = -p {- € ijdjV + di N L) , h m = -pN, 

with the convention for Levi-Civita tensor density tn = 1. Here, all the spatial indices are raised 
and lowered by 5^. A further note on this specific choice of decomposition is about the £ 2 /t 2 
coefficients: With this coefficients, at every step the flat space limit I — > oo, £/t — > 1 will be clear. 

Unlike the flat space case, <p is not gauge invariant anymore. In fact, under the gauge trans- 
formations Sftnv = V^C^ + ^uCfi where can be decomposed as = (Co, ~ e ijdjC + the 
components of h^ v transform as 



k N L = p ( « + Co + 7 « ) , <*<AT = 2^ ( Co + 7 Co 



Again, from the linearized Bianchi identity, V jU ^ I/ = 0, we know that there should be three 
independent gauge-invariant combinations constructed out of the (derivatives of) six scalar fields. 
By inspection, one can find these combinations, but the quickest way would be to look at the 
independent components of the gauge-invariant tensor . This led us to the following four 
gauge-invariant functions: 



f = 



2 lr 11 



. + X— t N 



V = 



-N 



V 2 N 



+ X ~ 2V 2 N L ~j(N- 2V 2 iV L + x ) + ^n\, a = - f (c - V 2 t?) , 



and a relation between them coming from the Bianchi identity 



f-p+- 



p — q = 0. 



(15) 



In terms of these, the components of the linearized Einstein tensor can be found as 



Goo = ~Y^ 2 ^ 
t 



— (dip + e ik d k a) 



(Sij + did^j q - didjp - (e ik d k dj + e jk d k d^j a 



The linearized curvature scalar follows as 



^ = i(^-v 2 / + p) = ^v 2 (/- P ), 



where in the second line we used the Bianchi identity. 
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Using the above, the Einstein-Hilbert action can be reduced to the following form: 



Ieh 



-a h^ v Ql = - I d 6 x 



As in the flat space case, computations get a lot simpler if the higher-derivative parts of the 
Lagrangian are organized in such a way that is replaced by some gauge-invariant combinations. 
This can be done again upon use of the self-adjointness of the involved operators as follows: 



L 2ct+(3 



(2a + /3) 



J d 3 x ^gh, u (>□ - VV + ^gr) Rl = J <Px ^gRl 



For the j3 term, one has 



P 



$ / d A x V-gKv aSf - -xgTRL 



4 1 d* x 



iPh^QZ-^KRL 



After organizing R^ u (|14[) into a form where the indices fi and v in the covariant derivatives stay 
at the far left, and using the Bianchi identity, V ' Jd^ = ^> one arr i ves at 



l 



d 3 x 



u ( -2^gf + \r\ + 1 vsf 



Note that, had we not done this and instead computed h flu OQ£ u directly, putting the result into 
an explicitly gauge-invariant form would be somewhat time-consuming. Not worrying about the 
correct canonical dimensions for the fields, one can collect all the parts computed above to end up 
with 



2 



2/3 

a + 



^fRL + - e (fV 2 f+p 2 + a< 



£ 3 

+ (2a + P) ^R 2 L 



a 2 + <tVV + p 2 + pV 2 p + (V 2 /) ' 



+ £ lR L V 2 f -^R L p-pV 2 f 



The flat space limit of this action gives ([6]). In this form, not all the fields are independent: After 
defining tp = V 2 /, and using the Bianchi identity (fT5j) . we can further simplify the action to 



I 



1 



2/3\ t 



i 2 J * 



(-tpp + p 2 ) + (2a + p)-(<p- V 2 p) ' 

+/3 ^3 ~ pV2p ~ tifV 2 p ~ tplp - (ppj | + I a , 



where the a field decouples from the rest 



(16) 



(17) 



For vanishing a and f3, cosmological Einstein theory does not have any propagating degrees of 
freedom just like its flat space partner. For generic a and p, there are 3 degrees of freedom. Recall 
that a minimally coupled scalar field with the correct canonical dimension is in the following form: 



I 



d A x yf^gida^d^^ + m 2 3> 2 



d 3 x 



f -^m 2 $ 2 
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Therefore, after rescaling a — > j^a in (|17p . one finds the mass of the a field as 



a 

1 



1 



12a 



Kp l 2 P 



(18) 



For generic a and /3, unlike the flat space case, diagonalizing the <p, p action is highly nontrivial. 
But, there are various ways to see the basic oscillators in this model. One such method is to Fourier 
transform the fields just in the x space and then consider the zero two-momentum limit. That 
would be equivalent to dropping the V 2 terms in the action. Note that this construction does not 
change the number of degrees of freedom, of course as long as V 2 (field) is not the lowest order 
term. Another way is to directly study the equations of motion. We shall employ both of these 
methods below. 

A. Masses from the nonrelativistic limit 



Apart from the decoupled a part, the generic a, j3 theory (|16j) reads in the nonrelativistic limit 



as 



/ - - / d 3 x 



a + 



20\ t 



\ t t f 

-tptp + p 2 ) + (2a + 0) -3 Cp 2 + p— [f -p> 2 - tptp - cpp 



To decouple the fields, first note that 2a + (3 = f + 8 "+ 3/3 , and rescale tp as p> — > jip to get the 
action 



/ - - / d 3 x 



2p\ t 



« + — 1 - (-PV + P 2 ) + j ^ ( <f 2 - ^ + 4p 2 - i,v 



(8a + 3/3) t 3 / 2 2,p> r 
+ 4 pi* 



Then, define a new field as $ = p — 2p, which leads to the decoupled actions for the $ and <p fields. 
As the spin-2 helicity partner of the a field, the $ action is exactly like the a action with the same 
mass m g (fTSjl : 



/* = ? / d 3 x 



e 3 e\p i 2] 



and the spin-0 mode has the action 



(8a + 3/3) 



d A x 



t 6 



1 t f 24a 6/3 



(8a + 3/3)£V I 2 £ 2 

which after putting into the canonical form by rescaling <p — > jzp> yields the mass 

2 1 4 / 3a + 



k (8a + 3/3) I 2 V8a + 3/3y ' 

In the 8a + 3/3 = case, the ip field freezes out and m 2 , matches the result of [l[ obtained with the 
help of an auxiliary field, not via canonical analysis. For generic a and /3, in accordance with the 
analysis of [l|, one can introduce two auxiliary fields to rewrite the action ([T]), but decoupling of 
the scalar mode from the spin-2 mode is not immediately clear. This is done in Appendix B. 
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B. Equations of motions in the BHT case 



The above nonrelativistic analysis reveals the canonical structure of the generic a, f3 theory. 
But here let us consider the relativistic equations of motion for the 8a + 3/3 = case. Dropping 
the a field in (1161). we have 



2 



+ t -p (p 2 - p^ 2 p -p 2 - t(pv 2 p - tpip - cppj I . 



It appears that there are 2 degrees of freedom in this action (which would conflict our earlier result, 
and the result of [1]), but this is a red herring, there is only a single degree of freedom. A quick 
way to see this is to look at the Hessian matrix, % — 



H 



4P 



dqidqj 

2 -2t\ 
-2t 4 ) 



t 



Since detH = 0, there is a constraint in the model. Therefore, "velocities" ip and p cannot be 
separately expressed in terms of the canonical momenta 



dC /3t 5 / . „ 2 2 • 



dC /3i 3 , n . 



One can use the Dirac's constraint analysis method to obtain the Hamiltonian for this singular 
Lagrangian, but here it suffices to consider just the field equations. Taking the variations with 
respect to (p and p yield 



Sip : 



m 2 „t 2 



t 6 / 

•'' p - -3 (2ip + tV 2 p + p 



W 9 ° 



t 5 (ip - V 2 p) - 2t 4 p 



0, 



and 



Sp : 



m 2 t 



(tp - 2p) 



2P 



V 2 [ip- V 2 p+^p + -p 



t 



t 3 {2p -tip-ip) 



0. 



By inspection, and with a hint from the field equations which give Rl = 0, one observes that 
<p = V 2 p and the other equation reduces to 



t 



1 



t 



<p--<P + V<p]-p[m a - ls ]<p = 0, 



1 



e- 



which is not yet in the canonical wave equation form in dS. To put in the canonical form, 



□ - m z 



0, rescale cp — > ip/t to obtain 

-ip + ~ip + V 2 p) - ^ m 2 g v = 0, (□ - m 2 ) p = 0, 



which is exactly like the a field. 
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IV. CONCLUSIONS 



We have studied the canonical structure of the linearized quadratic gravity models in an explic- 
itly gauge-invariant way for both flat and dS backgrounds in three dimensions. In flat spacetime, 
the general action is decoupled into three harmonic oscillators. After considering the signs and 
various limits of the parameters k, a, /3, the BHT case is singled out as the unique unitary and 
nontachyonic theory (namely, a regular massive free spin-2 field, not a higher-time derivative one), 
while the others are all higher-derivative Pais-Uhlenbeck oscillators. Sources are also added to the 
theory, and Newtonian potentials for both static and spinning particles are calculated. Moreover, 
we have computed the weak field limit of the circularly symmetric spacetime. We extended our 
flat space analysis to include the gravitational Chern-Simons term and investigated the oscillator 
structure for the BHT limit: We have seen that in this limit the oscillators decouple with different 
masses, violating parity as expected. In dS, we have also found the most general action in terms of 
three gauge-invariant functions constructed from the (derivatives of the) components of the metric 
perturbation and carried out the decoupling of the fields in the nonrelativistic limit at the level 
of the action and in a relativistic form at the level of the field equations. For future work, to 
go beyond the free field level and introduce nonlinear ities, such as O (/i 3 ) and interactions, our 
gauge-invariant actions will be of great use. Another interesting point about the models that we 
discussed here is that, especially in (anti)-de Sitter backgrounds, for certain tuned values of the 
parameters novel phenomena such as partial masslessness or chiral gravity arise. These topics will 
be addressed in a separate work. 
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Appendix A: Spinning masses 



It is also of some interest to understand how spinning point particles interact in the generic 
higher-derivative model. This can be done as follows: First, note that the energy-momentum tensor 
for a massive (m) spinning (j) pointlike source is 



Too = m5^ (r - fi) , T* = h^djS^ (f- n) , 



Tij = 0. 



For two such conserved sources scattering amplitude was computed in [5j as 



4A 



d 3 x -2r;„ 



pa z + -□ 

K 



1-1-1 r 1-1-1 



fin z + -□ 

K 



T-T' 



to + 3/3) □ - -□ 

K 



From the nonspinning case, the only added part will be 



1 



AT' /3D 2 + -D T 



-i 



niO 



f5m? g 



( 



r — nj 



□ □ - m 2 g 



WW (r-r 2 ). 



After carrying out the space integrations, it reads 



1 



47^ /3D 2 + -□ T 



-ii0 



3\32 

"27T/3 



Kq (m g |n - f 2 \) , 
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for r*i ^ r 2 . Then, the total Newtonian potential energy, U = j4/time, becomes 

U = ^~ ( m i m 2 + 4m^u 2 ) #o (ra fl In. - r 2 |) - — mim 2 if (m s \n - r 2 |) . 

Since ji and j 2 could be of any sign, the part coming from the spin-spin interaction can be repulsive 
or attractive. In the BHT limit the last term disappears. 



Appendix B: The a, j3 theory with auxiliary fields 



Consider the quadratic Lagrangian ([I]) in three dimensions. Using two auxiliary fields 4> and 



f n V , one can rewrite it as 



c = I 

K 



where m 2 



" k(8o+3/3) 



9 

and mk 



r - rc uv -w + ^ + r A [ rUv - f 



linearized 



hT + / 



2 4 

^. After linearization around flat spacetime, we have 
2 ,o 1 



re/3 



7 /xi^ 



k (8a + 3/3) 



4k/3 



hi/ 



r 



For 8a + 3/3 = 0, decouples, and f uu can be eliminated to yield the action describing spin-2 field 
with a Pauli-Fierz mass[l]. But, for generic a and /3, one has to find a way to decouple <fr , 
and h[i U keeping in mind that there should be a kinetic term for the <p field. This is possible by 
rescaling h nu , but we have not pursued this [20]. 



Appendix C: Linearized field equations in the de Sitter background 



In the body of the text, we worked mostly at the level of the action. To check our results at 
the level of the field equations, some of the computations in this Appendix are needed. The trace 
of the linearized field equation is 



ia + 3/3) OR L + 



6 (4a + /3) 



Rl = 0, 



where g^Q^v = was useci - Without further ado, let us list the results of somewhat tedious, 
yet relevant computations: 



°<?oo = — 



I 3 

2£ 3 



V 2 / + \ V 2 / - V 2 V 2 /) - \ V 2 p - | V 2 / - -— R L 



t 5 



= ^ [p+- t p- V 2 p - ±p - -V 2 f\ + ^jdj [&+-&- rv - ^o- \ . 



3 
— i 

t 



2 



^ didj U + ip + ^p- v 2 p - lv 2 p - |v 2 / 



2£ 3 



(tiAdj + £jk&kdi) (^ + j^ + j2 & - v2( 7 - 7 v2cr ) • 
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Qj^ u , and Rl, computed in the body of the text, together with the Bianchi identity (|15|) . and the 
above results are sufficient to study the field equations. 
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